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Ch ■ Abstract 

^ , 

^ I In this note we explicitly solve the Lindblad equation for a system 

of three spins with a three-body interaction, coupled to the environ- 
ment by bath operators that inject or absorb spin carriers. We exem- 
plify the properties of this solution in the context of a simple instance 
Q^ ' of Feynman's quantum computer in which a two-qubit program line 

is executed, applying the VNOT primitive to a one-qubit register. 
\0 . 

O ■ PACS: 03.67.Lx, 03.65.Yz 

o. 

1 Introduction 

Much effort has been spent in recent years on the problem of quantum trans- 
rS I port along a 1-dimensional spin chain with or without the presence of some 

C^ ■ interaction with an external environment [H [21 |3l HJ |5]. There are, however, 

only few examples of exactly solvable open quantum systems[6l [7]. In this 
note we present the explicit solution of the Lindblad equation for a simple 
(the simplest, indeed) spin chain in which the interaction between two neigh- 
bouring sites is mediated by another spin located on the link that connects 
them. Such a three-body interaction is familiar in the context of lattice 
gauge theories [8] but we prefer to examine its role in the timing mechanism 
of Feynman's model of a quantum computer |9]. 

In Feynman's machine, the application of computational primitives to the 
register is timed by a quantum clock: a single "spin up" travels along a spin 
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chain and is coupled to the register in such a way that, when it moves from 
site L to site L + 1, the L-th operation is performed on the register. In an 
ideahzed mode of operation, one supposes that the timing excitation travels 
ballistically along the chain. 

In this note we begin to explore the effect of excitations being injected/absorbed 
by the environment. Having in mind the Imry-Landauer (TUj model of con- 
ductance, we model the environment through the bath operators suggested 
in section III.B of reference 4. Guided by the idea that computation is a 
nonequilibrium capability of a nonstationary system, we pay more attention, 
in this note, to short-time transient behaviour than to time-asymptotic prop- 
erties. 

The paper is organized as follows. In section II we introduce the model and 
review its ballistic behaviour. In section III we study the coupling of the 
system to the environment via a Lindblad equation. Section IV is devoted to 
the discussion of examples. Conclusions and outlook are in section V. The 
spectrum of the Lindblad operator and a nonequilibrium stationary state 
(NESS) are studied in the appendix. 

2 The model: ballistic behaviour 



Figure 1: The complete system: the endpoints of a spin chain of length A^ 
interact with two reservoirs Bi, Br. The interaction between sites L and 
L + 1 is mediated by an ancillary spin a through the unitary operator U{a). 

Figure [1] represents the basic functional unit of Feynman's model of a 
quantum computer [9l ITTl [T2l [T3l fll] . The spin 1/2 systems located at the ver- 
tices L and L + 1 of a graph, that we indicate by r(L) = (ti(L), T2(L), T3(L)) 
and r(L + 1) = (ri(L + 1), r2(L + 1), r3(L + 1)), are part of the clocking 
mechanism or cursor. 

We will make the simplification of neglecting all other cursor sites (A^ = 2). 
This will be crucial in making our model explicitly solvable. The effect of a 



long clocking chain is, by the way, quite well understood [ISl E] and we wish 
to concentrate here on other effects. 

In an idealized mode of operation (to be reviewed below), the vertices L and 
L + 1 support an excitation that, when travelling along the oriented edge 
(L,L + 1), applies the unitary transformation U{a) to a register two-level 
system a = (o"i, o"2) cts). The inverse transformation f/^(cr) is similarly asso- 
ciated with the oriented edge (L + 1, L). 

We will suppose in this section that, setting, in our simple model, L = 1 and 
T±{j) = (ti(j) ± iT2{j)), for j = L,L + 1, the above degrees of freedom are 
coupled by a Hamiltonian of the form: 

H = -lf/(a) r+(L + l)r_(L) - ^U\a) r^L + l)r+(L). (1) 

We will take, by way of example, U{a) = e~*2'<^i, namely a rotation by an 
angle ^ around the axis 1. 

With this choice, the Hamiltonian H is best described in a basis of simul- 
taneous eigenvectors of (ai, T2,{L), T2,{L + 1)). We will list the vectors of this 
basis in lexicographic order (with "1" preceding "-1") as: 

|ei) = |ai = l,r3(L) = l,r3(L + l) = l) 
I 62 ) = I ai = 1, r3(L) = 1, r3(L + 1) = -1 ) 

I eg ) = I ai = -1, TsiL) = -1, r;{L + !) = -! ). 
In this representation the only non vanishing matrix elements of H are: 
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( 62 \H\ 63) = {cj \H\ eg ) = --e*2 = ( 63 |i/| 62 ) = ( ce \H\ 67 ). 

The matrix ||( Cj |iJ| ejt )|| j=i,...,8 is easily diagonalized and exponentiated 

leading to a quite explicit expression of the time evolution operator Z{t) = 
exp{—itH). In particular, for any initial density matrix pin we can ex- 
press the solution of the quantum Liouville equation dpt/dt = —i[H,pt] as 

Pt = Z{t)pinZ^{t). 

In this note, where not otherwise stated, we will always suppose that the 
state evolves from the following initial condition: 

Po = ^(/ - o,){l + n{L)){l - n{L + 1)). (2) 



This is the typical initial state of a Feynman machine: the excitation is 
initially located at the beginning of the spin chain, and the register a is in 
the state corresponding to the input (as = —1 in our specific example). 
Under this initial condition, the probability that T'i{L) = 1 is given by: 



ns{L,t) =Tr 



Pt- 



i + n{L) 



the probability that r3(L + 1) = 1 is given by: 

I + T3{L+l] 



n3(L + l,t) =Tr 



Pi- 



cas i- 



sm 



While the "spin up" excitation bounces back and forth between the two 
cursor sites L and L + 1, the register spin a is acted upon by the operators 
U and W , so that its Bloch vector evolves as: 

mi(t) = Tr[pt di] =0; 

m2(t) = Tr [pt 0-2] = sin (t/2) sin(^); 

ms{t) = Tr [pt (T3] = sin (t/2)' (1 - cos(?9)) - 1. 

We adopt, in most of the examples that follow, the value i) = 7r/2. We refer 
to the primitive U{tt/2) as to y/NOT, because by applying it twice the initial 
o"3 = — 1 state gets flipped into a^ = +1. 

For the sake of comparison with the behaviour studied in the following sec- 
tions, it is interesting to look at the Bloch diagram of the register spin, 
shown in figure |21^a). The motion of the system is periodic: at times t = 
{2k — l)7r, /c = 1, 2, . . ., the computation y/WOT is completed and the state 
of the register is on the Bloch surface, i.e. its the von Neumann entropy is 
zero. At times t = 2fc7r the computation turns out to be undone, the state 
is brought back to the surface of the Bloch sphere, the clocking excitation is 
brought back to its initial position. 

Of particular interest for the considerations that follow is the graph (figure 
[2]^c)) of the von Neumann entropy of the register: 



S^'^^t) 



l+r{t) ^^fl + r{t) 



1 



^In" 



r(t) 



where r{t) = ^m2{tY + ms{ty. Since the system is a closed bipartite quan- 
tum system, its overall entropy is zero, and the two subsystems, register and 
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S^(t) 
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Figure 2: Idealized mode of operation; i) = 7r/2. (a) A parametric plot of 
{m2it) , ms(t)) for < t < 27r shown inside the unit circle, (b) m2(t) (solid 
line) and m3{t) (dashed line) as functions of time, (c) The von Neumann 
entropy of the register as a function of time, (d) Solid line: ns{L,t), dashed 
line: n^lL + l,t). 



cursor, are isoentropic. The fact that there are times at which the intended 
computation (the rotation by 7r/2 of a spin initially pointing "down" ) is com- 
pleted with certainty is related to the possibility of bringing, at such times, 
S'^^^{t) to its initial value 0. Any source of noise, therefore, might blur such 
a certainty. 

As a final remark of this section we wish to give a hint to the development 
of the entanglement during the evolution under the Hamiltonian ([1]) of the 
initial condition (EJ. To this end we study the mean value of the chirality 
operator: 

X = ^-r(l)xr(2), 

whose interest as an entanglement witness is discussed in reference 16. It is 



easy to check that: 

Tr (x pt) = 2 cos(^/2) sm(t) 

so that it is only for -^ = that the state of the system attains the maximum 
value 2 of mean chirality compatible with bipartite entanglement. 

3 Lindblad equation 

The model outlined in the previous section is so simple that we can try to 
study the effect on it of interaction with the environment. 
Following references 4 and 5, we describe this interaction in the Markovian 
approximation and adopt a Lindblad equation approach, in which the fol- 
lowing bath operators act on the cursor spins: 

Li = Ve(l-/i)r+(L), L2 = v/e(l+/i)r_(L), (3) 

L3 = ^e{l + fi)T+{L + l), L4 = Ve(l-^)r_(L + l), 

where the coupling parameter e and the asymmetry parameter /i satisfy the 
conditions e > and — 1 < /i < 1. The asymmetry parameter /z models a 
possible difference of the chemical potentials of Bl and Br. We refer the 
reader to the second section of reference 17 for a discussion of the range of 
validity of this approach and of the possible degree of control of the baths 
and of the bath-system interaction. 
The density operator is supposed to evolve according to the equation: 



dpi 
dt 



1 ■* '^ 

-I [H, Pt] - 2 E {^1^^' ^*} + E ^^P^^r (4) 



i=i i=i 



The main aim of this note is to show that the Cauchy problem posed by (|2]), 
(|3]) and @ can be given an explicit solution, parametrized by {}, e, p. This 
is shown in full detail in the appendix. The full spectrum of the Lindblad 
matrix is, furthermore, provided there for the reader interested in initial con- 
ditions other than (jj]). 

In this section we study some physical properties of this solution. 
We observe, first of all, that we can give a quite explicit expression for 
the joint probability distribution of T^i^L) and r^i^L + 1); setting pt{j,k) = 
{ Gj \pt\ Cfc ) and using the solution (see the appendix) of System{+,+) we 
have: 



i. Probability that T3(L) = 1 Ar3(L + l) = 1 and probability that T^iL) = —1 
Ar3(L + l) = -l: 



Tr 





// + r3(L)\ / 


Pt 

Tr 


[ 2 j[ 


1 + 


r[ 2 

4e2 (1 - ^2) 



— e 



'l + Ts{L + l] 



I-TsjL + l] 

2 

1 + 4e2 (1 + /i) 



Pi(l,l)+Pi(5,5) 



-4te 



+ e 



4(l + 4e2) 4(l + 4e2) 

_2te^l^ (2e (1 + /x) cos(t) + sin(t)) 
l + 4e2 ■ 



We observe that these two probabilities vanish for e = (conservation of 
the number of clocking excitations in the Hamiltonian regime); for e > the 
bath operators can kill the single clocking excitation introduced in the initial 
condition (|2]) or introduce a spurious one. 

a. Probability that t^^L) = — 1 A r3(L + 1) 



1: 



Tr Pt 



I-Ts{L)\ fI + rs{L + l) 



Pt{3,3) + Pt{7,7) 



1 + Ae^l + ^if l + Ae^l + f,y 



-2te 



4 (1 + 4e2) 



4(l + 4e2) 



(l + 4e2(l + /i)')cos(t) 
2(1 + 4e2) ■ 



in. Probability that T^iL) = 1 A r^i^L + 1) = —1: 



Tr Pt 



I + n{L)\ fI-r;{L + l] 



p<(2,2)+pi(6,6) 



l+462(l-/i)'^^_4t.l + 46Ml+/i) 



+ e 



4(l + 4e2) 4(l + 4e2) 

24, (1 + 4e2 (1 - p^)) cos(t) - 4ep sin(t) 
2(l + 4e2) ■ 



The fact that the above joint distribution of t^{L) and t^{L + 1) does not 
depend on the parameter d justifies Feynman's statement that "zt turns out 
that the propagation of the cursor up and down this program line is exactly 
the same as it would he if the operator [U] were not in the Hamiltonian" [9]. 



This statement needs however a quahfication because of the ■(9- dependence 
of some components of the solution of System{+, +) and System^—, — ) ex- 
hibited in the appendix: there are, in fact, observables of the cursor (other 
than T-i{L) and tz{L + 1)) whose distribution depends on the primitive being 
apphed to the register. For instance, for the current operator, one has: 

Tr p,^(r_(LK(L + l)-r_(L + lK(L)) 

= '- (a(2, 3) - p,(3, 2) + p,(6, 7) - p,(7, 6)) 

cos I — 



2 J 4(l + 4e2) 
^\ 2e-2*^ (2e/i cos(t) + (1 + Ae^{l + y)) sin(t)) 



+ ^°^'.2; 4(l + 4e^) 

As to the study of the Bloch diagram of the register, we observe that, setting 
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m(t) = 2j]pi(j,j+4), 

it is: 

m2(t) = Tr [pt (T2] = Im {m{t)) , 
m^{t) = Tr [pt as] = Re (m(t)) . 

Under the initial condition ([2]) and under the assumption that the four eigen- 
values oji defined in the appendix are distinct, we can give an explicit expres- 
sion of the function m{t) in terms of the solution of System{+, —). 

Having the analytic solution of the master equation allows us to make some 
considerations about the asymptotic behaviour of the system as well. 
For example, the relaxation time of the register is determined by 

1 



max{Re (a;,) , i = 1, . . . , 4} 

The relaxation time of the cursor, on the other side, is (2e)~^, as it is 
possible to evince from the explicit expression of the matrix elements of 

System{+, +). 
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Other aspects of the system quantitatively accessible through the analytic 
solution are the asymptotic entropies of the overall system and of its compo- 
nents. It is easy to check that, for the register, it is: 

lim S""»(t) = log(2). 

For the overall system it is: 

lim Tr [-pt log {pt)] = log(8) - 4eV + O(e^). 

t—^oo 

Setting p1^^ = Tr„{pt), it is easy to check that the asymptotic entropy of the 
cursor 

lim Tr [-pr log [pD] = log(4) - 2eV2(cos(0) + 1) + O (e^) 

depends on the three parameters. 

As a final remark of this section, we point out that the mean chirality is 

given by: 

Tr{xpt) = 8 1m{pt{2,7)). 



4 Examples 

Figures El^d) and lU^d) show, for different values of the parameter /i, the 
marginal distributions of T3(L) and r^lL + 1) characterized, respectively, by 
the expectation values: 

Tr[p,(/ + r3(L))/2] = 

_ 1 + Ae'^{l- fi) _2t,(l + 4e^(l+yu))cos(t) -2e/isin(t) 
2(l + 4e2) +'' 2(l + 4e2) ' 



Tr[p,{I + niL + l))/2] = 

_ 1 + 4e2(l + /i) _2t^ (1 + 4e2(l + p)) cos(t) - 2ep sin(t) 
2 (1 + 4e2) ^ 2(l + 4e2) ' 

By comparison with figure |2t^d) it is immediate to notice, in both frames, the 
damping effect of the coupling parameter e and, in figure Ht^d), the polarizing 
effect of the asymmetry parameter p: the left reservoir injects spins up in the 
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Figure 3: e = 0.1; j2 = 0; ^ = 7r/2. (a) {m2{t) , m^it)) <t < QOtt; (b) as 
in figure 1. (c) Solid thin line: the entropy of the overall open system; solid 
thick line: the entropy of the cursor; dashed line: the entropy of the register, 
(d) Solid line: Tr [pt (1 + r3(L))/2]; dashed line: Tr [pt (1 + r^iL + l))/2]. 



chain while the right bath absorbs spins up, thus favoring a positive current. 

In this sense, when p y^ 0, the environment acts as a "battery" [TO] . 

A significant effect of the coupling with the environment can be seen in the 

entropy landscapes of figures [3]^ c) and|l](c). In both examples, the entropy 

of the overall system increases monotonically and the entropies of the two 

subsystems (register and cursor) differ form each other; the behaviour of the 

entropy of the cursor in the two cases is qualitatively similar. 

What is interesting, in the short-time transient we are mainly interested in, is 

how differently the entropy of the register behaves in the absence or presence 

of a chemical potential difference between the reservoirs. 

For p = (reservoirs at the same chemical potential) the entropy of the 

register increases almost monotonically, following the trend of the entropy of 

the cursor (figure [3]^c)). 

For p = —1 (figure 111(c)), around time n, the entropy of the register gets 
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Figure 4: e = 0.1; /i = —1; ?? = 7r/2. Frames (a),(b),(c) and (d) as in figure 

El 

close to the value 0, as it happened in the isolated system (figure EJ^c)). 
The presence of the "battery" is able to push, at t = vr, the Bloch diagram 
closer to the boundary of the Bloch sphere (figure Hl^c)), which, we remind, 
corresponds to the register being in the desired output state with certainty. 

As the coupling with the environment increases, the polarization of t{L) 
and t{L + 1) increases accordingly (figure [5]^d)). The presence of a "strong" 
battery results in the onset of a stationary current, witnessed by the uniform 
rotation of the register (see figure EJ^a)). The different relaxation times of 
the register and the cursor become evident, together with the dependence 
of the entropies of the cursor and of the overall system on the interaction 
parameters e, /i and ■&. 

As a final example of this section, we exhibit in figure 6 the evolution of the 
mean chirality Tr {x Pt) under Lindblad dynamics with the same parameters 
as in figure 4. 
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Figure 5: e = 1, /x = — l,-;? = vr/lO. (a) Parametric plot of {m2{t) , m^^t)) . 
(b) Re{ui) < Re{uj2) < Re^Us) < Re{cOi) as functions of {}. (c) To compare 
with [3]^ c) andm^c); this time the entropy of the register is larger than the 
entropy of the cursor, (d) To compare with[3l|^d) andlH^d). 
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Figure 6: e 
time. 




1; •& = 7r/2. The mean chirality as a function of 



5 Conclusions and outlook 

Our interest is in the effect of noise deteriorating the operation of the ma- 
chine. We model by differences of chemical potential the action of a battery 
trying to compensate for the noise. The examples discussed in the previous 
section suggest that, as far as the short-time transient is concerned, the dis- 
ruptive effects of the interaction with the environment can be mitigated by 
a proper tuning of the (macroscopic) reservoirs. 

The toy model we discuss in this paper is an oversimplified version of the 
systems discussed in references 2 and 3. With respect to the two references 
just quoted, however, our simple model adds a peculiar difficulty: the links of 
the chain do not carry numerical hopping parameters, but additional quan- 
tum degrees of freedom, those of the register on which the computation is 
performed. This added ingredient sets the Lindblad equation @ out of the 
family of problems explicitly solvable by third quantization [7]. 
In this note we used the entropy of the register to measure the "degree of 
certainty" (fidelity) with which the output state is reached; but figure \5\ sug- 
gests also to use the register as a probe to visualize what happens in the 
cursor when it interacts with the environment in more complicated ways. 
For example, one could investigate the onset of Anderson localization in a 
cursor with noisy hopping parameters. 

An explicit solution of the model of figure [H that is a system with an ar- 
bitrarily long cursor spin chain, would be most desirable. Our study of the 
structure of the master equation gives us some hint about how to tackle this 
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general problem: as far as the dynamics of System{+, +) and System{—, — ) 
is concerned, the solution scheme of third quantization |7j can be still applied. 
This means that observables of the cursor are already accessible. What re- 
quires a different approach is the solution of the subproblems System{+, — ) 
and System{—, +), which contain all the information about the register and 
the register/cursor interaction. 

Appendix: Spectrum of the Lindblad operator 

The analysis of the system (|4]) of 64 equations is made simple by the fact 
that the commutation relations 

[a^,H] = 0; [a,,L,] = 0, j = l,...,4 

and the adoption of the basis {\ e^), A; = l,...,8} make it possible to split 
the problem into 4 uncoupled systems of 16 equations. Setting, in what 
follows, pt{j,k) = ( Cj \pt\ Ck ), the four systems can be described in the 
following way: 

System{+, +) relates matrix elements of pj between vectors of the basis {|efe)} 
both belonging to the eigenvalue +1 of ai and is of the form: 



dpt{j,k) 



dt 

m— 1,... ,4 
n^l,...,4 



E ^;Jm,nP*K^), l<J,fc<4. 



The 16 X 16 matrix M++ turns out to have characteristic polynomial: 

det(M++ -xl)=x{l + {x + 2ef) (x + 4e) (x + 2e)^ 
x'^ + 4ex + 3e^ - e/i + - j (x^ + Aex + 3e^ + e/i + - 
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Under the initial condition (|2]) only the following six components of the so- 
lution of System{+, +) do not vanish identically: 

P^^'^'^= 8(l + 4e2) +^ 4(l + 4e2) + 

4(l + 4e2) 8(l + 4e2) ' 

1+46^(1 -/if ^,(._2e) 1 + 4^6/i + 46^ (1 - /i^) , 

^*^' ^" 8(l + 4e2) +^ 8(l + 4e2) + 



8(l + 4e2) 8(l + 4e2) 

8(l + 4e2)" 



p^(3, 3) = ^^Jf^Wf (1 - e*(^-2^) - e*(-'-2^) + e^^"^^)) ; 



a(4, 4) = ^ - (P<(1, 1) + Pt{2, 2) + p,(3, 3)) 



p<(2,3) = pi(3,2) = 

_ itf/2 /^ ^e^ , t(^-2.) ^-2e(l + /^) f(_^_2,) ^ + 2e(l+p) 
V2(l + 4e2)^ 8(2e-2) 8(2e + z) 

System{—, — ) relates matrix elements of pj between vectors of the basis {|efc)} 
both belonging to the eigenvalue —1 of ai; it is of the form: 

^^^ = E M-k-m,nPtim, n), 5 < J, fc < 8. 

TTl — 5,...,8 

n=5,...,8 

The 16 X 16 matrix M can be obtained from M++ by changing i) into —i) 
and has, therefore, the same eigenvalues. By the same argument it is easy 
to check that under (|2]) the only non vanishing components of the solution 
of System{—, — ) are: 

Pt(5,5) = pi(l,l), pi(6,6)=pt(2,2), 
pt(7,7) = pt(3,3), pi(8,8)=pt(4,4), 
pt(6, 7) = e-^Vt(2, 3), pt(7, 6) = e*V(3, 2). 
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System{+, — ) relates matrix elements ( e^- \pt\ e^ ), where | e^ ) belongs to the 
eigenvalues +1 of cxi and | e^ ) belongs to the eigenvalues —1 of ai. 

dptU, k) ^ ^ M^^^^^^^m, n),l<j<4,5<k<8. 



dt , , 

m — 1,...,4 
n^5,...,8 



1 , - ^4 



det(M"'" — xl) = — (x + 2e) (x — cui) {x — 102) {x — C03) {x — U4) 

[1 + 8 (2x^ + 16ex^ + (1 + 44e2) x^ + 4 (e + 12e^) x + e^ (4 + 18e^ - A^^)) + 
- Se^ ((1 + /i^) cos(^) + 2ifxsm{d))]^. 

Here (a;i,a;2,C(;3,a;4) is any ordering (for instance by increasing real part) of 
the four eigenvalues 

u±^± = — 2e± 

± ^■\/4e2 - 1 ± v^TTiee^ - 8eV + Se^ (cos (^) + /i^ cos(^) + 2i/isin(^)) 

V2 

of M"*"" that actually contribute to the evolution of the non vanishing com- 
ponents of the solution under the initial condition (jj]). 
In terms of the eigenvalues uj±^±, we can write: 

p,(l,5)=pi(4,8) = 

^ 6 A ^^ Jfi - l)e---^ + (p + 1) (-1 + 8e\fi - 1) + 4e(/i - 2)u, - 2u;J) 
4^ Y[i<k<4 {ujj - ujk) 



1 A ^.,, (2e + c.,) (1 + 4e^(p - 1)^ - 4e(p - 2)a;, + 2a;|) 
^*^ ' ^ 4^^ n^f^^^(^.-^fc) 



g-J" 1 dl 

p,(2,7) = -pi(3,6) = 

We notice, for completeness, that the spectrum of M~^~ includes also the 
following eigenvalues (which do not contribute to evolution under the initial 
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condition (^), each with multiphcity 2: 



C±,± = -2e ± - — A/e2i'? (4e2 - 1) ± 2e''^\J&'^t^ (/i - 1 + e*'' (/i + 1))'. 

System^—, +): relates matrix elements ( Cj \pt\ Ck), where | Cj ) belongs to the 
eigenvalues —1 of cxi and | e^ ) belongs to the eigenvalues +1 of cxi. It is of 
the form: 

m = 5,...,8 
n = l,...,4 

for 5 < J < 8, l<k<A. 

The 16 X 16 matrix M '" can be obtained from M^ by Hermitian conjuga- 
tion. The solution of System^—, +) is, obviously, obtained from the solution 
of System{+, — ) by Hermitian conjugation. 

We observe that, as only M+"'' and M have the eigenvalue 0, the initial 
condition po converges, as t — )■ oo, to a Nonequilibrium Stationary State Pness 
represented by the density matrix having the following non vanishing matrix 
elements: 



Pness(l, 1) = Pness{4:, 4) = Pness{5, 5) = Pness{8, 8) 

l + Ae^l-pf 



1 + 4e2 (1 - /i2 
8 (1 + 4e2) 



Pness (2, 2) = Pness (6, 6) 
Pnessy^i 'j) ^ Pnessyiy ') 



8(l + 4e2) 

l + 4e2(l + /if 

8 (1 + 4e2) 



Pnessl2,OJ — PnessK'J^ ^) — Pnessl' 5 6j — Pness\'^^ <) 



ie/ie*''/^ 



2(l + 4e2)' 
It is immediate to check that this state has entropy: 

/ 4,2 2 X 26ptanh-i(^^^^^) 
^— = log(8) - log 1 - -^^4^ , Vi+4. (i+My _ 
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